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POPULAR SUMMARY 

Clouds play an important role in our climate system. For example, clouds can 
reflect part of the solar radiation and reduce the energy arriving at the surface, thus 
cooling the climate. They can also absorb radiation emitted by the surface and reemit a 
portion back, thus warming the climate. It is critical to understand the properties of 
clouds in order to predict the climate. Satellites have become the most important tool in 
studying clouds. Most satellites use passive remote sensing techniques (which observe 
the Earth by receiving the radiation reflected or emitted by it). However, in recent years, 
active remote sensing satellites (which observe the Earth by sending out radiation 
actively and then observing the reflected signal) have become an essential component of 
the Earth observing system. 

The Geoscience Laser Altimeter System (GLAS) is an active remote sensing 
system. As the first satellite lidar mission with global coverage capability, GLAS was 
launched in January 2003 on board the Ice, Cloud and land Elevation Satellite (ICESat). 
Since its launch, GLAS has contributed significantly to the study of cloud properties. For 
example, using GLAS data, we can obtain the exact cloud top height and accurately 
distinguish cloudy from clear areas in the polar regions, which are very difficult to 
accomplish using passive remote sensing techniques. However, GLAS has limitations 



too. Laser beams emitted from GLAS can only penetrate clouds to a limit of a few optical 
depths (optical depth is a measure of transparency. The larger the optical depth, the less 
transparent it becomes). Hence, only the optical depths of thin clouds can be retrieved by 
GLAS. Our study showed that about two-thirds of the clouds were too thick to be 
penetrated by the laser. 

Is it possible to obtain thick cloud optical depths from GLAS and other space borne 
lidar instruments? Our work proves the answer is yes. We did this by using the sunlight 
that is reflected by clouds and which then reach the GLAS sensor. The reflected sunlight 
registers as a noise to the GLAS system because it coincides with the reflected laser 
signal. For GLAS to retrieve cloud properties, this noise has to be subtracted. However, 
the reflected sunlight carries information about clouds too and this is why passive remote 
sensing can retrieve cloud information. To use this information, we first calibrated the 
solar radiation received by the GLAS detectors. In so doing, we can treat GLAS as a 
passive remote sensing device and thick cloud optical depth can be obtained. Our method 
was tested against the thick cloud optical depth datasets retrieved from Geostationary 
Operational Environmental Satellite (GOES). The comparison showed that the average 
difference between the two datasets is 24%. As a result of, this study, thick cloud optical 
depths will be provided as a supplementary product to the existing GLAS cloud products 
in future GLAS data releases. For the first time, we will be able to retrieve passive and 
active remote sensing results simultaneously from one system. 



Retrievals of Thick Cloud Optical Depth from the Geoscience Laser 
Altimeter System (GLAS) by Calibration of Solar Background Signal 


Yuekui Yang 1 , Alexander Marshak 2 , J. Christine Chiu 3 , Warren J. Wiscombe 2 , Stephen 
P. Palm 4 , Anthony B. Davis 5 , Douglas A. Spangenberg 4 , Louis Nguyen 6 , James D. 

Spinhime 2 , and Patrick Minnis 6 


Journal of the Atmospheric Sciences 


'Goddard Earth Science and Technology Center, UMBC, Baltimore, MD. 
2 Goddard Space Flight Center, Greenbelt, MD. 

3 Joint Center for Earth Systems Technology, UMBC, Baltimore, MD. 
4 Science Systems and Applications Inc., Lanham, MD 
5 Los Alamos National Laboratory, Los Alamos, NM. 

^ASA Langley Research Center, Hampton, VA. 


1 



Abstract 


Laser beams emitted from the Geoscience Laser Altimeter System (GLAS), as well 
as other space-bome laser instruments, can only penetrate clouds to a limit of a few 
optical depths. As a result, only optical depths of thinner clouds (< about 3 for GLAS) 
are retrieved from the reflected lidar signal. This paper presents a comprehensive study 
of possible retrievals of optical depth of thick clouds using solar background light and 
treating GLAS as a solar radiometer. To do so we first calibrate the reflected solar 
radiation received by the photon-counting detectors of GLAS’ 532 nm channel, which is 
the primary channel for atmospheric products. The solar background radiation is 
regarded as a noise to be subtracted in the retrieval process of the lidar products. 

However, once calibrated, it becomes a signal that can be used in studying the properties 
of optically thick clouds. In this paper, three calibration methods are presented: (1) 
calibration with coincident airborne and GLAS observations; (2) calibration with 
coincident Geostationary Operational Environmental Satellite (GOES) and GLAS 
observations of deep convective clouds; (3) calibration from the first principles using 
optical depth of thin water clouds over ocean retrieved by GLAS active remote sensing. 
Results from the three methods agree well with each other. Cloud optical depth (COD) is 
retrieved from the calibrated solar background signal using a one-channel retrieval. 
Comparison with COD retrieved from GOES during GLAS overpasses shows that the 
average difference between the two retrievals is 24%. As an example, the COD values 
retrieved from GLAS solar background are illustrated for a marine stratocumulus cloud 
field that is too thick to be penetrated by the GLAS laser. Based on this study, optical 
depths for thick clouds will be provided as a supplementary product to the existing 
operational GLAS cloud products in future GLAS data releases. 
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1. Introduction 


The Geoscience Laser Altimeter System (GLAS) was launched on board the Ice, 
Cloud and land Elevation Satellite (ICESat) in January 2003 as part of the NASA Earth 
Observing System project (Spinhime et al., 2005a). GLAS observes Earth in two 
wavelengths, the 532 nm channel, which uses photon-counting detectors, and the 1064 
nm channel that uses analog detection. More sensitive to atmospheric signals, the 532 
nm channel is used as the primary channel for atmospheric products (Palm et al., 2002). 
Since its launch, GLAS has been providing data that contribute significantly to studying 
cloud and aerosol properties (e.g., Hart et al., 2005; Hlavka et al., 2005; Spinhirne, et al., 
2005b). However, the retrieved optical depths are limited to the relatively thin clouds 
that can be penetrated by the laser beam (< about 3). 

Prior to the lidar retrieval process, the reflected solar energy has to be subtracted as 
noise from the signals received by the photon detectors. However, Platt et al. (1998, 
2006) suggested that, if calibrated, the solar background can be viewed as a signal and 
used to retrieve cloud optical depths of dense clouds, thus completing the cloud probing 
capability of active remote sensing with lidar. The reflected solar energy is recorded by 
GLAS in units of photon counts. Calibration is needed to convert photon counts into 
radiances. One path to calibration is from the instrumental parameters that are measured 
in the lab. This method suffers from the uncertainties stemming from degradation or 
change of the instrument during its deployment. Valencia et al. (2004) proposed a 
method using collocated NASA Aerosol Robotic Network (AERONET) sunphotometers 
in calibrating the solar background of the ground-based micropulse lidars (MPL). 
Applying this method, Chiu et al. (2007) demonstrated encouraging results in retrieving 
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cloud optical depths for thick clouds. Their validations against other instruments show 
that retrieved cloud optical depths agree within 10-15% for overcast stratus and broken 
clouds. 

In this paper, we present a comprehensive study of three possible ways of 
conducting on-orbit calibration of the reflected solar radiation received by the photon- 
counting detectors of the GLAS 532 nm channel. Section 2 gives basic information on 
GLAS solar background signals. The three calibration methods are introduced in Section 
3. Section 4 demonstrates the validation of the calibration by comparing Geostationary 
Operational Environmental Satellite (GOES) and GLAS retrievals of Cloud Optical 
Depth (COD). A case study is presented in Section 5 to illustrate how bonus information 
can be obtained from calibrated solar background signal in addition to the results of 
GLAS active remote sensing. Our conclusions are stated in Section 6. 

2. GLAS Solar Background Signal 

The solar background component of a GLAS profile is derived by averaging the 
signal in the background portion of the profile where no laser signal is present. The 
background region used for GLAS signal acquisition is both prior to the laser beam 
reaching the atmosphere (about 80 km altitude) and well after the beam strikes the Earth: 
the equivalent of a range of approximately 40 km below ground (Palm, et al., 2002). The 
signal received by the GLAS photon counting detectors within these regions is 
overwhelmingly dominated by reflected sunlight (detector dark current is negligible).. In 
this study, we use the solar background signal determined with the background region 
above the atmosphere. 
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The GLAS solar background data are archived in the products GLA02 and GLA07, 
uncalibrated and calibrated lidar signal profiles, with horizontal resolutions of 40 Hz (175 
m) and 5 Hz (1.4 km) averaged over eight returns. These data are in units of raw photon 
counts and stored before the detector dead-time correction is conducted. Dead-time is a 
span of time immediately following the receipt of a photon during which the photon 
counting detector is unable to record the arrival of additional photons (Campbell et al., 
2002). The dead-time correction is performed by using a look-up table which contains a 
dead-time corrected value for each possible output from the photon counting channel 
(Palm et al., 2002). The dead-time corrected solar background photon counts are then 
used in the calibration process. 

Proportional to the corrected photon counts n (photon counts/bin) registered at the 
detectors, the radiance L (W-m'^sf'-pm' 1 ) of the reflected solar energy that reaches the 
GLAS instrument can be written as, 

L = Cn (1) 

where C is the calibration coefficient in the unit of W-m' 2 -sr' 1 -pm' 1 /(photon counts/bin). 
The calibration process is the practice of determining the calibration coefficient. 

The GLAS data used in this study come from the first campaign with full on-orbit 
operation of the instrument (termed L2A) that began on 25 September and lasted until 19 
November, 2003(Spinhime et al., 2005a). During this campaign, over two thirds of the 
clouds observed by GLAS were not penetrated by the laser. Hence, calibrated solar 
background will provide important complementary information to the GLAS active 
remote sensing products. In addition, some GLAS active remote sensing data collected 
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from later campaigns are with degraded quality due to technical problems with the lasers 
(Spinhime et al., 2005b), but the photon detectors that receive solar background signals 
have remained stable, therefore consistent COD retrievals can be expected from the 
properly calibrated solar radiation. 


3. Calibration Methods 

To reduce uncertainties in the calibration, it is best to employ multiple independent 
methods. Three methods are used in this study: (1) calibration with collocated 
MODIS/ASTER Airborne Simulator (MASTER) and GLAS observations; (2) calibration 
with collocated observations of deep-convection clouds by GOES and GLAS 
observations; (3) calibration from first principles using optical depth of thin water clouds 
over ocean retrieved by GLAS active remote sensing. 

3.1. Calibration with collocated GLAS and airborne observations 

The airborne observation data employed in this study are from the GLAS 
Validation Experiment executed with the high-altitude NASA ER-2 aircraft from NASA 
Dryden Flight Research Center in Edwards, California, in October 2003 (Hlavka et al., 
2005). Four instruments participated in the campaign: the MODIS/ASTER Airborne 
Simulator, the Cloud Physics Lidar (CPL), the Video Imaging System (VIS) and the 
MODIS (Moderate Resolution Imaging Spectroradiometer) Airborne Simulator (MAS). 
However, MAS was on the ER-2 for only 3 of the 7 GLAS missions and none of them 
were during daytime. So radiance observations from MASTER, which was on-board the 
ER-2 for all 7 missions, are used in this study. 
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MASTER was developed to support scientific studies by the Advanced Spaceborne 
Thermal Emission and Reflection Radiometer (ASTER) and Moderate Resolution 
Imaging Spectroradiometer (MODIS) projects (Hook et al., 2001). In sunlit regions, the 
radiance observed by ASTER’s 538 nm channel is close enough to the 532 nm channel of 
GLAS to be directly used in the calibration of the GLAS solar background. However, 
most of the October 2003 campaign flights were conducted at night. Among the 7 flight 
missions, only the one that took place in the early morning of October 24 th had sufficient 
sunlight along the track suitable for calibrating GLAS solar background radiances. 

Figure 1 gives the radiance image of the ASTER 538 nm channel observation (Fig. la) 
and the corresponding GLAS 532 nm channel attenuated backscatter image (Fig. lb) for 
the track at the time of the GLAS overpass. A special feature of the GLAS satellite is the 
ability to accurately point the lidar to within 50m of ground locations. Thus comparison 
of the satellite and aircraft data is possible. 

The size of the MASTER image (Fig. la) is 289x36 km 2 with a pixel resolution of 
50 m. The horizontal resolution of the GLAS image (Fig. lb) is 175 m. For calibration 
purposes, observations from both instruments need to be collocated both in space and in 
time. Collocation in space is done with the nearest neighbor technique and the accuracy 
is within 50 m crosstrack. Then three MASTER pixels along the track are averaged to 
match the size of the GLAS sampling distance. Due to the speed differences between the 
two platforms, most of the pixels collocated in space are not collocated in time. The 
image scan time for the MASTER was 22 minutes and 24 seconds, while for GLAS it 
was 42 seconds. To minimize the ensuing uncertainties, we limit the pixels used in this 
study to those within 5 minutes time differences between the two observations. Figure 2 
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gives the flowchart of the calibration process. 

Figure 3 shows the calibration result. The two pixel clusters in Fig. 3 represent the 
collocated clear and the cloudy pixels. As is customary in instrument calibration (e.g., 
Nguyen et al., 2001), the regression line is forced through the origin. Physically, this is 
because corresponding to zero photon counts, the solar background radiance must be zero 
as well. Equation (2) gives the regression, 

L = 6.62 n (2) 

So the calibration coefficient derived with this method is C-6.62 W-m‘ -sr' -pm" /(photon 
counts/bin) and the one-sigma error of the slope is 0.06. 

When the radiances corresponding to the solar background photon counts of each 
selected point are calculated using Eq. (2) and then compared to the MASTER 
observations, the mean difference is 4.1% with a standard deviation of 3.3%. A variety 
of sources may contribute to this calibration uncertainty, including the remaining space 
and time difference and the size difference between the pixels. 

3.2 Calibration with collocated deep-convection clouds observed from GLAS and 
GOES 

The second approach employs collocated GLAS and GOES visible channel deep 
convection observations. Because of the single line nature of GLAS images, exact 
collocations in time are rare between the spatially collocated GOES and GLAS pixels. 
However, because of the horizontal homogeneity of deep convective systems and 
relatively weak fluctuation in radiances reflected from very thick clouds, we can assume 
that small differences in time and space between the selected GLAS and GOES pixels 
would not cause significant bias in the calibration results. After collocated deep 
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convective observations are selected, the solar energy differences between the GOES 
visible channel (0.65 pm) and the GLAS green channel (0.532 [im), and the view angle 
differences between the observations from the two instruments must be taken into 
account in the calibration process. 

The following criteria for selecting deep convection points are employed: (1) GOES 
10.7 pm channel brightness temperature < 205K and its 3x3 pixel standard deviation < 
IK; (2) the 3x3 pixel standard deviation of the GOES 0.65 pm channel raw count < 3% 
of the central pixel raw count; (3) GLAS reported cloud top height >10 km; (4) 
temperature at cloud top from GLAS products < 208K. In addition, to have sufficient 
sunlight, the solar elevation angles for the selected points had to be > 12°. Twenty-one 
collocated deep convection points are found with GOES- 10 and GOES- 12 data during the 
GLAS L2A campaign period. Table 1 lists these points; entries having longitudes < - 
105° are from GOES-10. As shown in the table, due to the differences in the observation 
strategies of the two instruments, the selected pixels are collocated in space, while 
differences exist in the observation time. 


Figure 4 gives the flowchart of the calibration process with this method. Since the 
pixel sizes of GLAS (175 m) and GOES (4 km) are different, the GLAS solar background 
signal is averaged to match the GOES pixel resolution. The GOES-10 and 12 radiances 
are first calibrated to the corresponding Terra MODIS 0.63-pm channel using the 
methods described by Minnis et al. (2002). They are then adjusted to a 532 nm 
wavelength from the original 630 nm measurements with the following equation: 


L =L x Msn 

^532 ^650 * 




630 


( 3 ) 
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9 1 

where M 5 3 2 =1869 W-m' -pm' is the solar spectrum irradiance at 532 nm and M 6 3o=1641 
W-m" 2 -pni 1 is the solar spectrum irradiance at 630 nm (ASTM, 2000). 


As listed in Table 1, the view angles of the selected GOES observations are always 
away from nadir. These view angle differences are taken into account by using the 
Angular Distribution Models (ADM) developed by the Clouds and the Earth's Radiant 
Energy System (CERES) Inversion Working Group (Loeb et al., 2001) with the 
following equation: 


L 0 (d 0 ) - Le(9 0 ,(p 0 ) x 


ADMM) 

ADM g (G 0 ,cp 0 ) 


(4) 


where Zo(0o) is the radiance at nadir for solar zenith angle 0o, £e(0o,<Po) is the radiance at 
view angle 0 for solar zenith angle 0o and solar azimuth angle cpo, ADMo(Qo) is the ADM 
value at nadir for solar zenith angle 0o, and v4DMe(0o,(po) is the ADM value at view angle 
0 for solar zenith angle 0o and solar azimuth angle qpo. 

Figure 5 shows the calibration results with this method. As with the first method, 
the regression line is forced through the origin, 

L = 6.36 n (5) 


The calibration coefficient derived with this method is C=6.36 W-m‘ 2 -sr' 1 -pm' 1 /(photon 
counts/bin) and the one-sigma error of the slope is 1.63. 


A variety of factors can cause uncertainties to the calibration coefficient determined 
with this method. First, the adjustment process of the GOES’ radiance data may bring 
inaccuracy. For example, the CERES ADMs are derived with broadband observations. 
Even though the bandwidth of GOES’ visible channel is fairly wide (0.52 - 0.72 pm), 
remapping GOES’ off-nadir observations to nadir with CERES ADMs can still result in 
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biases. Second, the collocation process can be another error source. For example, time 
differences exist in the collocated pixels. The largest difference (that corresponds to the 
smallest, <150 W-m^-sr' 1 , radiance in Fig. 5) is for Point 17 in Table 1, which is for a 
0.24-hour (14.4 minutes) time difference. Even though the reflected solar energy for 
deep convective clouds is usually stable, the differences in observation time can still 
cause uncertainty to the calibration coefficient. 

3.3. Calibration from first principles 

The third approach takes advantage of the active remote sensing results from 
GLAS. It involves three steps: (/) determine the reflected solar radiances using radiative 
transfer calculations for GLAS-retrieved thin cloud optical depths as input; (if) select 
points with the lowest solar background at each cloud optical depth; (mi) derive the 
calibration equation from the calculated radiances and the measured solar background 
photon counts. The thin cloud optical depth is a standard GLAS product (GLA1 1) that is 
retrieved from analysis of the lidar backscattered signal (Spinhime et al., 2005b). By a 
thin cloud, we refer to a cloud that does not completely attenuate the lidar signal 
(generally, with optical depth < about 3). 

Figure 6a plots the solar background photon counts versus the GLAS GLA1 1 cloud 
optical depths over ocean for the GLAS L2A campaign. Due to the uncertainties in the 
phase function of ice clouds, only single layer liquid clouds (with cloud top height < 

3500 m) have been selected. The plot contains around 18,000 points that have solar 
zenith angles ranging from 60° to 70°. As shown in the figure, there is a wide spread in 
the values of solar background photon counts that correspond to each retrieved cloud 
optical depth. A variety of reasons, such as surface reflectance variability, differences in 
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aerosol loading, cloud microphysics and uncertainties in the retrieval process, can result 
in different solar reflectance for clouds with the same optical depth. It would not be 
practical to determine the at-sensor solar radiances for all the points. However, the lower 
boundary of the scattered points, which is marked as a thick line in Fig. 6a, represents the 
observations with the lowest solar background that corresponds to the lowest surface 
reflection and the least aerosol loading. It is feasible to calculate the reflected solar 
radiances corresponding to these observations using radiative transfer models. Figure 6b 
shows GLAS retrieved cloud optical depths versus the solar background photon counts 
for these points, which are selected through the following procedure: (/) the observations 
are binned with an optical depth interval of 0.2; (//') points with optical depths smaller 
than 0.02 and larger than 0.8 are excluded to keep only the most reliable GLAS 
retrievals; and (///) four points with the lowest solar background photon counts in each 
bin are selected. We limited the data to only warm water clouds to avoid additional 
uncertainties related to the scattering phase function. 

Simple ID radiative transfer calculations are conducted with the DISORT (Discrete 
Ordinates Radiative Transfer Program for a Multi-Layered Plane-Parallel Medium) 
model (Stamnes et al., 1988) under the following assumptions: (/) since the points are for 
the lowest surface reflectance, the wind speed according to the Cox and Munk (1954) 
model is assumed to be small (5 m/s); (//) the aerosol optical depth is assumed to be 0 
(lowest aerosol loading); and (3) the cloud effective radius i? e ff is assumed to be 10 p,m. 
The uncertainty caused by this assumption is studied and presented in Fig. 8. 

Figure 7 gives the flowchart of the calibration process with this method and Fig. 8 
shows the calibration results. If the regression is forced through the origin then 
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L = 6.35 n. 


( 6 ) 


2 

Hence, the calibration coefficient derived with this method is C=6.35 W-m' -sr' 
pm' 1 /(photon counts/bin) and the one-sigma error of the slope is 0.84. As mentioned 
above, the radiances are calculated by assuming a cloud droplet size of 10 pm. The 
vertical error bars in Fig. 8 give the uncertainties caused by a typical range of droplet 
effective radius (6 pm < i? e ff <16 pm). As seen in the figure, the uncertainties are small 
(with a maximum value of 4.3%) and do not affect the calibration coefficient 
significantly. 

3.4. Combination of the three methods 

The calibration coefficients (6.62, 6.36 and 6.35) derived from the three methods 
agree well with each other, within 4.2%. To finalize the results, all the points used in the 
three methods are consolidated together and plotted in Fig. 9. 

Equations (7) and (8) give the linear regression relations between the solar 
background x and the reflected solar radiance L with and without forcing the regression 
line through the origin: 


L = 6.38 n 

(7) 

L = 6.34« + 0.84 

(8) 


Based on this, we determine the final calibration coefficient C=6. 38 W-m’^sr’^pm' 
V(photon counts/bin). The one-sigma error of the least-squares slope is 0.05. 


4. Comparison of COD Retrievals from GLAS Solar Background and from GOES 
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4.1. Data 


Once calibrated, the reflected solar background signal received by the GLAS 
photon counters can be employed in retrieving the optical depths of thick clouds. The 
retrieval process is straightforward. Since the retrieval is based on the information from a 
single channel, we have to make an assumption about the value of the effective radius of 
the cloud droplets. In this study, we use i? e ff = 10 pm as a baseline value. As will be 
shown later, the possible bias caused by this assumption is usually within 10%. The 
lookup tables for the retrieval are built with DISORT (Stamnes et al., 1988) calculations 
with a 2° solar zenith angle interval and a 0.1 cloud optical depth interval. 

To test the validity of the COD retrievals from GLAS solar background, we 
compare them with the independent retrievals from GOES. The collocated GLAS and 
GOES observations of deep convective clouds used in the calibration process (Section 
3.2) are excluded from comparison. The cloud properties from the GOES data are 
determined with the Visible IR Solar-IR Split Window Technique (VISST) (Minnis et al., 
1995, 1998), which categorizes clouds into water, ice, and super cooled liquid water 
phases. To simplify the comparison, we only use clouds with water or super cooled 
liquid water phases over ocean. 

Given the large region covered by GOES, a significant amount of spatially 
collocated points between GLAS and GOES can be found. However, the time 
differences between the two observations could be large. For this study, we use only the 
spatially-matched GLAS and GOES data points that occur within 15 minutes of each 
other. All together, 741 points were found that satisfy the above-mentioned requirement. 


14 



Figure 10 gives the distribution of time differences between the two observations for the 
selected points. 

Another problem in comparing GOES and GLAS COD retrievals is the different 
spatial resolution. The GOES cloud optical depth is taken from an approximate 16x16 
km 2 area centered on the GLAS point, while the GLAS footprint is 175 m. Consequently, 
the GLAS solar background signal has to be averaged over 92 points to ensure maximal 
overlap between the two retrievals. Since the area used to obtain the mean GLAS data 
values (175 x 16,100 m 2 ) has a different spatial size and shape compared to the GOES 
retrieval footprint, significant discrepancies are to be expected between the two retrievals 
for inhomogeneous cloud fields. 

As a first approximation, for the GLAS solar background retrievals we assume that 
the surface is not reflective. For solar zenith angles within the range between 50° and 
80°, the Cox-Munk model tells us that the ocean nadir reflectance is about 0.5-2% 
depending on the wind speed. We assume this reflection is insignificant for comparison 
with the GOES retrievals, especially for optically thick clouds. 

4.2. Results of comparison 

Figure 1 1 shows the results of comparing the cloud optical depths retrieved from 
the GLAS 532 nm solar background and from GOES. As seen from Fig. 1 la, there is a 
wide scatter of points with a bias towards higher COD retrieved from GOES. As 
mentioned above, the main source of discrepancy here is the difference in spatial 
resolution of the two datasets; this is especially true for highly inhomogeneous clouds. 

To illustrate, we calculated the standard deviation of GOES cloud optical depth within an 
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approximate 24 km x 24 km area and the standard deviation of the 92 GL AS points 
corresponding to the GOES retrieval footprint. These standard deviations represent the 
amount of cloud horizontal inhomogeneity. The smaller the standard deviations are, the 
more likely that the observed clouds with the two instruments share the same properties 
and the closer the retrievals should be to each other. Indeed, if the standard deviations 
are limited to 25% of their corresponding mean values, a much better correlation with 
essentially no bias between the two retrievals is achieved (Fig. 1 lb). On average, the 
relative difference, which is the mean absolute difference between GLAS and GOES 
COD over the mean of GOES COD, is 24%. And the relative root mean square (RMS) 
difference, which is the RMS of the difference between GLAS and GOES COD over the 
RMS of GOES COD, is 28%. These differences are comparable to the 25% RMS 
difference between MODIS optical depth retrievals using the VIS ST and optical depth 
retrievals based on measurements taken at the surface (Dong et al., 2008). 

Figure 11c illustrates examples of the radiance distributions for the two points 
highlighted in Fig. 11a, one with large standard deviation in GLAS data (labeled “L” in 
Fig. 11a), and the other with small standard deviation (labeled “S” in Fig. 11a). As seen 
from panel c, a small shift in space for the large standard deviation case would result in a 
large difference in radiance and hence a large difference in the retrieved COD, while the 
small standard deviation case does not have this problem and hence provides a good 
match between the GLAS and GOES retrievals. 

Two additional factors of the GLAS COD retrieval process may affect the 
correlation between the retrievals from the two instruments. The first factor is that the 
one-channel retrievals from GLAS solar background assume a fixed cloud droplet 
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effective radius (here 10 pm). The uncertainty arising from this assumption is illustrated 
in Fig. 12a. The upper and lower bounds of the retrieved COD are determined by 
assuming an effective radius of 6 pm and 16 pm. Based on the retrieved GOES effective 
radii, the 6-16 pm range covers 82% of the data. As shown in Fig. 12a, the RMS 
retrieval errors resulting from the uncertainty in effective radius is 7%. 

The other factor that can cause errors in the GLAS retrievals is the uncertainty in 
the calibration coefficient of GLAS solar background. As discussed in Section 3, the 
difference among the calibration coefficients derived from the three individual methods is 
within 4.2%. Here we assumed a 5% uncertainty in the calibration coefficient, which led 
to the error bars plotted in Fig. 12b. The error values are generally larger than those 
caused by the effective radius uncertainty with the root mean square of 15%. Obviously, 
the larger errors are for optically thicker clouds. 

To better understand the total error resulting from the uncertainties in both effective 
radius, R e ff, and calibration coefficient, C, we assume that both uncertainties are normally 
distributed (see insets in Fig. 13a and b), with a mean of 10 pm and a standard deviation 
of 3 pm for and a mean of 6.38 W-m' 2 -sr'' -pm' '/(photon counts/bin) and a standard 
deviation of 2.5% for C. (Note that while simulating the sensitivity to effective radius, 
we used a truncated normal distribution rejecting R e ff below 6 pm and above 16 pm.) 

The distribution of the retrieved COD values is calculated using a straightforward Monte 
Carlo procedure picking randomly realizations of i? e ff and C. Figure 13 shows two 
examples for thicker (Fig. 13a) and thinner clouds (Fig. 13b). For the thicker cloud, the 
solar background was 32.1 (photon counts/bin) while for the thinner cloud it was 18.7 
(photon counts/bin). The calibration coefficient, C=6.38 W-m' -sr' -pm' /(photon 
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counts/bin), and effective radius, R t ff=\0 pm, lead to a COD=37 and 1 1, for the thicker 
and thinner clouds, respectively. With normally distributed uncertainties in i? e ff and C, 
the resulting COD distribution has a mean of 37 and a standard deviation of 4 for the 
thicker cloud, and a mean of 1 1 and a standard deviation of 0.6 for the thinner one. This 
translates to 1 1% and 6% one-standard-deviation-errors for the thicker and thinner 
clouds, respectively. 


5. A Marine Stratocumulus Case Study 

To further illustrate how GLAS passive remote sensing complements GLAS active 
remote sensing, we show an example that involves a thick marine stratocumulus cloud. 
The marine stratocumulus scene (Fig. 14a) was observed by GLAS on 1 November 2003. 
The scene, which extended over 900 km, was taken over the southern Pacific Ocean from 
35.13°-43.29°S and 84.30°-85.80°W. The cloud deck is optically thick and the standard 
GLAS active remote sensing was unable to retrieve its optical depth. However, this 
information can be obtained using solar background signal. Figure 14b shows the 
retrieved COD field. The average COD for this scene is 1 1, which is typical for marine 
stratocumulus clouds. 

With the empirical Eq. (11) from Minnis et. al. (1992) derived for marine 
stratocumulus, its geometrical thickness, Ah, can be estimated based on COD, x, as 

Ah = 0.452x^. (9) 

The average geometrical thickness of the clouds in this scene is thus ~260 m. The cloud 
top height is determined by GLAS active remote sensing and is a standard GLAS product 
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(GLA09). In addition to the cloud top height, Fig. 14c also shows the cloud base height 
determined by subtracting cloud thickness A h from cloud top height. As a result, if the 
empirical relationships between A h and x are reliable, passive remote sensing can 
complement the active remote sensing by determining cloud base heights for clouds too 
thick to be penetrated by laser beams. 

6. Conclusions 

It has been suggested by Platt et al. (1998, 2006) that solar background count rates 
in spacebome lidar returns can be used to infer cloud optical depth as long as they are 
properly calibrated. In this paper, we examined three possible ways of calibrating the 
reflected solar radiation that reaches the GLAS’s 532 nm channel photon-counting 
detectors. In so doing, we turn solar background radiation, which has been so far 
regarded as noise to be subtracted in the retrieval process of the lidar products, into a 
signal that could be used in retrieving the optical depth of optically thick clouds, which 
can not be penetrated by the GLAS lasers. Three independent calibration methods are 
investigated: (a) calibration with collocated airborne and GLAS observations; (b) 
calibration with collocated deep-convection clouds from GOES and GLAS observations; 
(c) calibration from first principles using optical depth of thin water clouds over ocean 
retrieved by GLAS active remote sensing. The main results are the following: 

(1) The calibration results from the three methods agree well with each other and 
the differences among the calibration coefficients are within 4.2%. Consolidating all the 
data used in the calibration, we determined the final calibration coefficient to be 6.38 
W-m^-sr' 1 -[Am' 1 -(photon counts)" 1 . 
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(2) Cloud optical depths retrieved from calibrated GLAS 532 nm solar background 
radiances match those from the GOES satellites fairly well when the inhomogeneity of 
the cloud field is accounted for. The retrievals from the two instruments have a 
correlation coefficient 0.87 with essentially no bias. On average, the difference between 
COD retrieved from GOES and GLAS overpasses is 24%, a value similar to the 
difference between optical depths derived from satellite and surface instrumentation. The 
GLAS root mean square retrieval errors resulting from effective radius uncertainty are 
about 7% while the errors from possible calibration uncertainty are on the order of 15%. 

(3) The retrievals have been demonstrated for a GLAS scene with marine 
stratocumulus clouds too thick for the GLAS laser to penetrate. In addition to cloud top 
height retrieved from GLAS active remote sensing, we used the GLAS calibrated solar 
background signal to retrieve cloud optical depth. As an example, we then converted 
cloud optical depth into cloud geometrical thickness using an empirical relationship 
derived for marine stratocumulus (Minnis et al., 1992). This, combined with the direct 
lidar measurement of cloud top, allowed us to estimate cloud base. 

Based on this study, optical depths for thick clouds will be provided as a 
supplementary product to the existing operational GLAS cloud products in future GLAS 
data releases. The methods presented in this paper, even though implemented for GLAS, 
can be used to calibrate solar background signals for other space borne lidar instruments, 
such as the Lidar In-Space Technology Experiment (LITE) on space shuttle Discovery 
and the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard 
CALIPSO. We understand that CALIPSO, as a part of A-train, has MODIS onboard 
Aqua flying only 15 sec apart. Furthermore, CALIPSO itself has the Wide Field Camera 
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(FWC) that takes measurements at 645 nm, and is designed to match the Aqua MODIS 
instrument's channel 1 . However, for current and future missions without the advantages 
that CALIPSO has (e.g., ICESat II), the methods studied in this paper provide examples 
to follow. 
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Fig. 1. (a) The MASTER 538 nm image from the October 24, 2003 flight off the west 
coast of California. The thick white line in the middle of the image represents the GLAS 
track. The double arrow line marks the region that has time differences less than five 
minutes between GLAS and MASTER observations, (b) The corresponding GLAS 532 
nm attenuated backscatter image. 
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Fig. 3. Calibration of GLAS solar background signal with collocated MASTER 
observations. The calibration equation and the uncertainty in the slope are shown in the 
upper-left comer of the figure. Total number of points is 450. 
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Fig. 4. Flowchart of calibrating GLAS solar background signal with the collocated GLAS 
and GOES deep convection observations. 
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Fig. 5. Calibration of GLAS solar background signal with the collocated GOES deep- 
convection cloud observations. 
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Fig. 6 (a) GLAS retrieved cloud optical depth vs. solar background photon counts for 
water clouds. The thick black line (approximately) represents the lower boundary of the 
scattered points, (b) The selected points used in the calibration process. 
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Fig. 7. Flowchart of calibrating GLAS solar background signal with the cloud optical 
depths retrieved from GLAS active remote sensing. 
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Fig. 8. Calibrating GLAS solar background signal with the thin cloud optical depths 
retrieved from GLAS active remote sensing (data product GLA1 1). The selected points 
correspond to the lowest values of solar background for each optical depth (see text for 
details). The calibration coefficient is derived by assuming a cloud droplet size with 
iJefplO pm. The error bars give the uncertainties caused by the range of a possible i? e ff (6 
pm</? e ff< 16 pm). 
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Fig. 11. Comparison of cloud optical depth (COD) retrievals from GLAS 532 nm solar 
background and from the GOES satellites. The GOES retrieval footprint is 16x16 km 2 
and the GLAS retrieval is derived from the mean radiance of 92 175 m resolution data 
values, (a) Results for all 741 available GOES points. Only retrievals with COD<100 
are plotted leaving 17 points outside the plot area. The two points marked in the figure 
represent cases with large (marked “Z”) and small (marked “S”) standard deviation of the 
GLAS data, (b) Results for the selected 73 points that have a standard deviation smaller 
than 25% of their mean values. Regression equations and the correlation coefficients are 
shown, (c) Radiance distributions of the GLAS data used in the calculations for the two 
selected points marked in (a). 
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(am to 16 (am) and (b) errors resulted from the uncertainty in calibration coefficient (5%). 






(a) 



(b) 

Fig. 13. Examples of the error analysis in the retrieved cloud optical depth (COD) for 
thicker (a) and thinner clouds (b). Insets show the assumed uncertainties in effective 
radius, R^n, and the calibrated background radiance resulted from the uncertainties in the 
calibration coefficient, C. R e s and C are assumed to be normally distributed. (Note that 
the small (R e ff< 6 pm) and large (i? e ff> 16 pm) values have been rejected.) Mean R e tn\0 
pm with standard deviation (std) of 3 pm, and mean C=6.38 W-m^-sr' 1 -pirf'^photon 
counts/bin) with std of 2.5% lead to mean COD=37 with std 4 and to COD=l 1 with std 
0.6, for the thicker and thinner clouds, respectively. 
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Fig. 14. A marine stratocumulus scene over the southern Pacific Ocean observed on 
November 1 st , 2003, (a) GLAS 532 nm backscattering image; (b) COD retrieved from 
GLAS 532 nm solar background at resolution 0.2 sec (1.4 km)); and (c) cloud top 
observed by GLAS mapped to the same resolution as in panel (b) and cloud base derived 
from an empirical equation (Minnis et. al. 1992). 
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Table 1. Collocated GOES and GLAS pixels during the L2A campaign 


Point 

Date of 
2003 

Latitude 

Longitude 

GOES 

GLAS 

Time 

(GMT) 

View 

Angle 

Time 

(GMT) 

View 

Angle 

1 

5-Oct 

9.53 

-71.46 

12.89 

11.95 

12.95 

Nadir 

2 

5 -Oct 

9.28 

-71.50 

12.89 

11.67 

12.95 

Nadir 

3 

10-Oct 

-4.06 

-63.76 

12.44 

14.60 

12.21 

Nadir 

4 

11 -Oct 

-11.20 

-67.66 

12.46 

16.04 

12.41 

Nadir 

5 

11 -Oct 

-11.70 

-67.73 

12.46 

16.47 

12.41 

Nadir 

6 

11 -Oct 

-11.95 

-67.76 

12.47 

16.69 

12.41 

Nadir 

7 

11 -Oct 

-12.20 

-67.80 

12.47 

16.92 

12.41 

Nadir 

8 

11-Oct 

9.09 

-161.65 

18.64 

32.75 

18.76 

Nadir 

9 

14-Oct 

6.80 

-73.96 

12.9 

8.20 

12.82 

Nadir 

10 

18-Oct 

17.08 

-84.21 

13.36 

22.42 

13.42 

Nadir 

11 

21 -Oct 

-11.20 

-72.70 

12.46 

13.58 

12.43 

Nadir 

12 

21 -Oct 

9.33 

-142.46 

17.14 

14.03 

17.17 

Nadir 

13 

21-Oct 

9.08 

-142.50 

17.14 

13.83 

17.17 

Nadir 

14 

21-Oct 

8.84 

-142.54 

17.14 

13.64 

17.17 

Nadir 

15 

21-Oct 

7.84 

-142.68 

17.15 

12.89 

17.17 

Nadir 

! 16 

21-Oct 

6.03 

-167.12 

18.65 

37.96 

18.79 

Nadir 

17 

29-Oct 

12.84 

-92.79 

13.38 

25.56 

13.62 

Nadir 

18 

4-Nov 

-9.21 

-40.79 

9.95 

40.97 

9.87 

Nadir 

19 

4-Nov 

-9.46 

-40.83 

9.96 

41.00 

9.87 

Nadir 

20 

4-Nov 

-9.71 

-40.86 

9.96 

41.03 

9.87 

Nadir 

21 

7-Nov 

9.10 

-71.25 

11.93 

11.57 

11.88 

Nadir 
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